Interference Effects in the 
Conductance of Multi-Level Quantum Dots 



^1" 
O 
O 

(N 

5h . 

< 

0\ 



C.A. Biisser, G.B. Martins, K.A. Al-Hassanieh, Adriana Moreo, and Elbio Dagotto 1 

'National High Magnetic Field Lab and Department of Physics, Florida State University, Tallahassee, FL 32306 

(Dated: February 2, 2008) 

Using exact-diagonalization techniques supplemented by a Dyson equation embedding procedure, the trans- 
port properties of multilevel quantum dots are investigated in the Kondo regime. The conductance can be 
decomposed into the contributions of each level. It is shown that these channels can carry a different phase, and 
destructive interference processes are observed when the phase difference between them is ±7r. This effect is 
very different from those observed in bulk metals with magnetic impurities, where the phase differences play 
no significant role. The effect is also different from other recent studies of interference processes in dots, as 
discussed in the text. In particular, no external magnetic field is here introduced, and the hopping amplitudes 
dot-leads for all levels are the same. However, conductance cancellations induced by interactions are still ob- 
served. Another interesting effect reported here is the formation of localized states that do not participate in 
the transport. When one of these states crosses the Fermi level, the electronic occupation of the quantum dot 
changes, modifying the many-body physics of the system and indirectly affecting the transport properties. Novel 
discontinuities between two finite conductance values can occur as the gate voltage is varied, as discussed here. 
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I. INTRODUCTION 

Among the several interesting potential applications of 
nanostructures and nanodevices is their role as small-size test 
ground for the analysis of the influence of many-body inter- 
actions in the physics of strongly correlated electronic sys- 
tems. These correlation effects are believed to be of crucial 
important in many bulk materials, such as high-temperature 
superconductors and manganites, but there is no universal 
agreement on what formalism is the most appropriate for their 
study. Verifying the accuracy of simple model Hamiltonians 
in nanoscale systems can provide an important test of the over- 
all validity of the theoretical formalism. Moreover, several 
compounds are believed to spontaneously form a complex pat- 
tern of nanoislands 1 and, as a consequence, accurate studies of 
a small number of atoms can produce interesting information 
even for bulk materials. 

In addition, new and novel effects can and do emerge from 
the quantum behavior of electrons in small environments, 
which are not obviously present in bulk systems. If quantum 
dots are considered as parts of electronic circuits, it is impor- 
tant to understand the new properties that may arise due to 
the discrete nature of the quantized states, the influence of in- 
teractions and spin, and the coherence of the wave functions 
in small systems that requires a fully quantum description of 
the problem. It is only recently that this complicated physics 
is starting to be unveiled and novel effects have already been 
found. A typical case is provided by the Kondo effect 2 , which 
in the bulk produces an increase in the resistance of a sample 
below the Kondo temperature 7k. On the contrary, in meso- 
scopic systems such as quantum dots, the same effect actually 
leads to a decrease in the resistance since the Kondo reso- 
nance provides a new channel for transport 3 ' 4 . 

The analysis of dots is not a mere academic problem. 
For example, transport properties at low temperatures across 
quantums dots (QDs) have received considerable attention in 
the last few years mainly due the possibility to use these struc- 



tures in quantum computing and spintronics. More specifi- 
cally, two coupled quantum dots (2QD) have been proposed 
to represent tunable qubit circuits 5 , where the spin interaction 
between the electrons in the dots defines the two states of the 
qubit. For all these reasons, the study of interacting electrons 
in nanosystems and its associated quantum transport are very 
important for progress in condensed matter physics in general. 

In several previous studies it has been usually assumed 
that, e.g., the 2QD have only one localized spin per dot, since 
only one level is assumed active in each dot. For the success of 
these proposals, it is important to understand the role played 
by other levels in the QDs - i.e. considering multilevel dots 
- providing part of the motivation for the work presented here. 
In addition, to properly perform these subtle studies it is cru- 
cial to carry out the calculations using unbiased techniques 
that do not assume properties of familiar states that arise from 
previous experience with bulk materials. Concepts such as 
Fermi liquids are and must be challenged in this context, since 
there is no reason to expect that nanoscopic entities will admit 
a similar description. Ideally, the calculations should make the 
less possible number of assumptions, and in this framework 
it is important to develop and use suitable numerical tech- 
niques to handle these complicated fully-quantum problems 
accurately. For this reason, the calculations discussed below 
were carried out with a powerful numerical method that has 
been recently developed. 

Other interesting investigations in this context have pro- 
posed the use of quantum dots as interferometers in the 
Aharanov-Bohm (AB) geometry with external magnetic fields 
applied 6 ' 7 . In general, previous efforts have considered a mag- 
netic flux enclosed in two possibles paths for the electrons. 
In these studies, the tunneling matrix-elements are defined 
with an explicit phase factor that represents the magnetic flux. 
These investigations have shown that narrow dips in the con- 
ductance could appear as a consequence of phase interference 
in the AB circuit, for electrons moving through different path- 
ways. Other investigations have also found conductance can- 
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cellations for the case of two quantum dots, individually con- 
nected to leads and among themselves 8 . This effect is present 
even using one-body interactions in the formalism, and it is 
caused by interference effects between two different paths in 
the geometry of the problem. 

One additional important motivation for the analysis dis- 
cussed in the present paper is to find alternative sources of 
interference among wave functions that can also lead to dips 
in conductance, in the absence of external fields and also in 
the absence of obviously distinct paths. Effects of this va- 
riety were recently addressed in multilevel dots by Kim and 
Hershfield 9 , where the consequences of having tunneling ma- 
trix elements with different signs in a multilevel QD and in an 
AB circuit were discussed. These authors showed that the 
matrix-element signs are important to determine the global 
phase of the electrons when they cross the QD. Also, arrays 
of odd-number interacting single-level quantum dots have in- 
teresting zeros in the conductance as the gate voltage varies 10 . 

In the investigations presented here, we consider a system 
composed of one or more quantum dots - each with two active 
levels and including the Hubbard many-body interactions - 
coupled to ideal single-channel leads. In gate-voltage regimes 
where there is only one electron in the QD, the usual Kondo 
physics 3 can be found in the low temperature regime. The 
coupling between the spin in the QD and the spins in the leads 
creates a resonance at the Fermi level that contributes to the 
transport. When a second electron enters in the QD a more 
interesting physics is observed since now the total spin of the 
system (singlet or triplet) plays an important role 111213 . One 
of the main results of our investigations is that the two levels 
of the dots studied here can act as different channels that indi- 
vidually contribute to the overall conductance with their own 
phase. This phase plays a role similar to that of a magnetic 
field in the AB interferometer, and it can induce constructive 
or destructive interferences. Below, it is explicitly shown that 
a phase diference of tt can exist between the two channels, 
even if there is no explicit difference between their hopping 
amplitudes and Coulombic interactions. Only a tiny differ- 
ence in the energies of the two levels is needed for the effect 
to develop, as well as a nonzero Hubbard coupling which is 
required to produce the Kondo physics. It should be remarked 
that this phase difference is different from the phase messured 
in some QD experiments 14 since it is internal to the QD. A 
magnetic field is not necessary to observe this interference. 

Others recent investigations have also focussed on sys- 
tems that present interference processes even without mag- 
netic fields or explicitly different paths. For instance, Hof- 
stetter and Schoeller analyzed lateral quantum dots near the 
singlet-triplet degeneracy point that occurs from the competi- 
tion of the Hund rule's coupling, that favors the triplet, and the 
nonzero energy difference between levels, that favors a singlet 
stated. A dip in the conductance was observed in this region. 
The study discussed in the present paper includes this singlet- 
triplet transition regime as a special case, but mainly focusses 
on a different variety of interference effects which are present 
even for a very small energy difference between the levels in 
the quantum dots. Our main picture is based on the notion that 
charge transport can occur with different phases when elec- 
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FIG. 1 : Schematic representation of the main system studied in this 
paper, consisting of one quantum-dot with two levels a and /3. The 
two possible paths carrying different phases & a and $,3 are shown. 

trons cross different levels of a quantum dot, rather than a 
singlet-triplet competition. Alternatively, our results can also 
be visualized as interference between two Kondo states: one 
with S — 1 and the other with S = 1/2. 

The organization of the paper is the following: In SeclITIthe 
discussion of the formalism, mainly the model and technique, 
is presented. In Sec.|lll|the case of one dot with two levels is 
studied in detail, showing the novel interferences in transport. 
The analysis of the integer-spin Kondo effect, that appears 
here in the interesting region of parameters, is in Sec. lIVI It is 
well-known that this type of states could enhance the conduc- 
tance near singlet-triplet transition s 11 ! 12 . A possible intuitive 
explanation of the effect is given in Sec.|V] An analysis of 
other cases where similar effects were found is in Sec. IVII 
and the conclusions are provided in Sec. IVHI Other effects 
analyzed in the paper include the possibility of observing lo- 
calized states near the dots. One such localized state in the QD 
can cross the Fermi level of the system without contributing 
to the conductance, but affecting the transport of other levels 
through the Coulomb interaction. This is shown to lead to 
discontinuities - as opposed to cancellations - in the conduc- 
tance vs. gate voltage, a concept that was not discussed before 
in the literature to our knowledge. 

n. ONE DOT WITH TWO LEVELS: FORMALISM 

The main system that will be analyzed in this paper, and 
especially in this section, consists of one quantum dot with 
several interacting levels, connected to ideal leads. It will be 
argued that - mainly due to the many-body interaction - in 
a transport process these levels can carry a different phase, 
as shown schematically in Fig. \l\ giving rise to a complex 
conductance pattern including destructive interferences (i.e. 
an "internal" Bohm-Aharonov effect appears to be dynami- 
cally generated). To model this system, the impurity Ander- 
son Hamiltonian with two levels (a and 0) is used. The por- 
tion corresponding to the isolated dot is 

Hdot = 51 U/2 n X(T n\ s + U'/2 2J n aa np a > + 

A— a,/3 a xjxj' 

— J y j C aa ' Cp a , Cfi a + 

crer' 
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where the first term represents the usual Coulomb repulsion 
between electrons in the same level, which for simplicity is 
considered equal for both. The second term represents the 
Coulomb repulsion between electrons in different levels (the 
U' notation is borrowed from standard many-orbital studies 
in atomic physics problems). The third term represents the 
Hund coupling ( J > 0) that favors the alignment of spins, 
and the last term is simply the energy of the states regulated 
by the gate voltage V g . Note that the level a is assumed to 
be separated from j3 by AV, and by modifying this param- 
eter an interpolation between one- and two-level physics can 
be obtained. The dot is connected to the leads (represented 
by semi-infinite chains) by a hopping term with amplitude *', 
while t is the hopping amplitude in the leads. More specifi- 
cally, 

-Hleads = 1 X! C K+1<7 + C ri(T c ri+la + h.C.] , (2) 

ia 

Him = t' ^ [ C \a ( C iOtr + CrOff) + h.c.] , (3) 

where c~l ia (c^ irT ) creates electrons at site i with spin a in the 
left (right) contact. Site "0" is the first site at the left and right 
of the dot, for each half-chain. The total Hamiltonian is 
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Note that for V g =-U/2 -U'+ J/2 - AV/2, the Hamiltonian 
is particle-hole symmetric. 

Using the Keldysh formalism 1 ^, the conductance through 
this system can be written as 



& -\t 2 G lr {E F )\ 2 \p{E F )] 2 , 



(5) 



where Gi r (E F ) is the Green function that corresponds to mov- 
ing an electron from the left lead to the right one, E F is the 
Fermi energy, and p(E F ) is the density of states of the leads 
(assumed the same left and right, for simplicity): 



-1 joj + if])- ^(uj + ir]) 2 -4t 2 
P(«) = ^ 



(6) 



To use Eq|5] |Gi r | 2 must be calculated. This quantity 
has two contributions, one corresponding to transport through 
path "a" and the other through path "/3". From the equations 
of motionii, which correspond to an exact expansion of the 
Green function™, it can be shown that 



G, 



gi t' G ar + gi *' G 



0r, 



(7) 



where gi are the Green function of the first left-contact and 
G ar (Gj3 r ) is the dressed Green functions to move from the 
state "a" ("/3") to the right-contact. 

The two terms in EqQ involving G ar and Gp r , can be eval- 
uated independently. To calculate the phase difference the fol- 
lowing expression was used 



i A$ = log 



G a r \Gj3r\ 



Gflr \G a 



(8) 



The conductance is proportional to |G/ r | 2 and, using Eqs. 
|5]Eland|S] it can be written as 
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(9) 



where a a and ap are the partial conductances given by 



o a = e -[t 2 git' p{E F )} 2 \Gar{E F )\ 2 , (10) 



[t 2 mt' p{E F )] 2 \G P r{E F ) 



(11) 



An exact cancellation of the conductance can occur when 
<j a = &f3 an d A$=±7r as a consequence of a destructive inter- 
ference process. In the same way, the situation with A<£>=0 
and <Ja=<jp can be considered as a constructive interference. 

The zero-temperature Green functions for this problem 
were calculated using an exact-diagonalization (Lanczos) 
method 19 to obtain the ground state and the Green functions of 
a small cluster containing the many-body interactions. This is 
followed by an "embedding" procedure to reestablish the rest 
of the leads 20,21 . The first step in this technique is to separate 
the system in two portions. One of them is a cluster containing 
the levels a and (3 (interacting region) and also including the 
first few sites left and right for each contact. The total num- 
ber of levels in the cluster is denoted by L. For one two-level 
quantum dot, this corresponds to L — 1 sites, since one of them 
has 2 active levels. Since Kondo physics will be important in 
this analysis, treating exactly not only the dot but also some 
sites belonging to the leads allows for a proper quantitative 
consideration of the spin-singlet Kondo "cloud". The other 
portion in the problem corresponds to the rest of the contacts. 
Denoting by g the exactly calculated Green functions (at T=0) 
of the cluster, the rest of the system can be incorporated using 
the Dyson equation G=g + g t G, where t are the matrix ele- 
ments that connect the cluster with the contacts, and G is the 
dressed Green function. 

To take into account possible charge fluctuations inside 
the cluster, a linear combination of g's calculated with dif- 
ferent number of particles is considered. For this purpose, 
g p =(l — p)g n + V 9n+i is defined, where g m is the Green 
function for m particles. With the dressed Green function G p , 
the total charge inside the cluster after the embedding is cal- 
culated using Q p =— 1/tt J2i l m ^ii( UJ )^ U! - O n tne other 
hand, the charge in the cluster before the embedding process 
was q p =(l — p)n + p(n + 1). These two quantities must be 
equal (Q p = q p ), defining a self-consistent equation for p. 

One of the important scales in this problem is the broaden- 
ing of levels of the dot due to the interaction with the leads. 
This energy scale, T=Trt' 2 p(E F ) 2 , defines the Kondo tem- 
perature for the case of a spin- 1/2 inside the quantum dot 2 . 
The coupling values used in this paper are (unless otherwise 
stated) U=0M, U'=2/3U, J=U - V, and t' = 0.2*, that 
allows for a fast convergence with the cluster size 21 . Fixing 
the values for the Coulombic interaction couplings is only for 
simplicity, to avoid exploring a vast parameter space. Consid- 
ering U 1 to be smaller than U is reasonable since the latter is 
expected to be the largest Coulombic interaction. The equa- 
tion for J is borrowed from atomic physics problems? 2 *, and 
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FIG. 2: Conductance for five different values of AV. The conduc- 
tance (in units of e 2 /h) is shown in solid lines and the phase differ- 
ence (divided by n) in dotted lines. When AV is the largest energy 
(AV=2t (e)) two structures similar to those found in a single-level 
QD, separated by AV, appear in the conductance. The case (d) is 
similar, but since AV is not large enough, an extra thin peak ap- 
pears around V g = — 1 (this peak is shown with more detail in Fig. 

When AV decreases, these structures began to interact giving 
rise to dips in the conductance, as shown in (b) and (c) for the cases 
AV=0.2t and AV=0.5t (AV=0.2t is studied in more detail in Fig. 
|4j- Finally, in the limit A V=0, the conductance dips are transformed 
into discontinuities. 



it gives a reasonably small value for the Hund interaction in 
units of U. In addition, we selected the ratio T/U to be close 
to the experimental values^. 



III. CONDUCTANCE CANCELLATIONS IN 
MULTILEVEL DOTS 

Using the Lanczos method to evaluate the Green functions 
of a small cluster followed by the embedding process de- 
scribed in the previous section, we have calculated the con- 
ductance, the charge of the QD, and the phase difference Eq|8] 
as a function of V g . Unless otherwise stated, we used the 
smallest possible cluster size L. For the case analyzed in this 
section this smallest size is L = 4 that correspond to two 
"lead" sites and two levels in the dot. In the last part of this 
section we analyze the convergence of the results increasing 
L. In Fig. |2 we show the conductance for several values of 
the energy separation AV between the two levels of the dot. 
At large AV (Fig. |2je)) the levels are well separated and the 
interaction between them is weak. The conductance shows 
just two wide peaks associated with the standard Kondo ef- 
fect for each level. Note that the width of these peaks is U, 
and with increasing temperature this broad peak will trans- 
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FIG. 3: Conductance (in units of e 2 /h) in a narrow range of V a and 
AV. The discontinuity at AV=0 is transformed into a dip for finite 
values of AV. 



form into two sharper peaks corresponding to the Coulomb 
blockade. This physics has been discussed at length in previ- 
ous publications and it basically amounts to the behav- 
ior of single-level QDs, which is natural at large AV. 

Far more interesting results are obtained by reducing AV. 
For instance, in Fig. |2jd) a thin central-peak is already ob- 
served. This peak seems to be the consequence of a construc- 
tive interference between the two possible channels, in a re- 
gion where a state with total spin <Sd=1 is formed, as will be 
discussed more extensively below. Upon further decreasing 
AV (Figs. |2jb) and (c) ) the two largest structures began to 
interact with the thin central peak, which increases its width, 
and a more complex pattern emerges. With the levels of the 
dot getting closer it is more likely to have them populated at 
the same time, and due to the Hund coupling the state with 
5d=1 has a high probability for the case of two electrons in 
the dot, giving rise to a spin-1 Kondo effect within the cen- 
tral peak. The other two peaks are associated with St>=1/2 
and 1 and 3 electrons in the dot. For AV<0.5t, dips in the 
conductance are clearly observed, and later we will show that 
these cancellations are related to interference processes. Fi- 
nally, at AV=0 (Fig. 0a)) several discontinuities appear in 
the conductance. It will be argued that these discontinuities 
are a consequence of the existence of localized states at the 
dot. These states are not of direct relevance for the transport 
of charge, but when the gate voltage moves them below the 
Fermi level of the leads they simply increase the charge of the 
quantum dot, which indirectly affects transport through the 
Coulombic interaction. Note that, for A V=0, the phase differ- 
ence between levels a and (3 is zero, since the Green functions 
G a 2 and G/32 are equal. Thus, no destructive interferences are 
observed, but discontinuities can still occur. 

In Fig. |5J it is shown how one of the conductance disconti- 
nuities is rapidly transformed into a dip as A V increases from 
zero. In Figs. |4j [5] and [6] we can observe in detail the case 
with AV=0.2t. Figure|4]shows the conductance together with 
the phase difference and the partial conductance. Clearly, the 
zeros of the conductance are in regions where A$=7r. Then, 
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FIG. 4: Conductance (in units of e 2 / h) (a), phase difference (divided 
by 7r) (b), and partial conductances (c) for AV=0.2t. The dips, indi- 
cated by "a" and "b" in (a), are in regions where the phase difference 
jumps to tv producing an exact cancellation when a a =<jp. 
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FIG. 5: Conductance (in units of e 2 /h), phase difference, and par- 
tial conductances for AV=0.2t for the "b" dip of Fig. |4j a )- Thicks 
arrows shows the gate potential where the conductance cancells, 
A<3> = 7r and o a — up. 



FIG. 6: Conductance (in units of e 2 /h), total charge of the QD, and 
charge at each level for the case AV=0.2t. The charge of the level 
a is indicated with a continuos lines, while the charge of /3 is shown 
with a dashed one. The central peak emerges in the regime where 
there are two electrons in the dot, in different levels. 
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FIG. 7: Conductance (in units of e 2 /h) and charges for the case 
AV=t as a function of the gate potential. These results should be 
compared against those of Fig. [6] 



from Eq|9] the dips must occur when the partial conductance 
a a and <jp take the same values. More detail for dip "b" is 
shown in Fig. [5] where thick arrows point to the zero in the 
conductance and the location where a a — erg. The dot charge 
for the same Ay is shown in Fig. [6] The central peak in the 
conductance occurs in the regime where two electrons are in 
the dot, in different levels. It can also be observed that the total 
charge at the QD n aa + np a ) does not show the char- 
acteristic sharp plateauxs of the Coulomb blockade regime: 
here the electrons are entering in the QD in a smoother way 
since the Coulomb interaction is comparable to the hopping 
amplitude t (Fig.|6jb) ). 

Similar results but for the case of AV=t are shown in Fig. 
[TJ In this situation, the large Ay well separates the levels 
a and [3, and as a consequence one of the levels is charged 
almost completely before some charge starts to enter into the 
other level (Fig. EI C ))- The latter occurs at the V g where the 



conductance peak is observed. 

To better understand the physics associated with the thin 
central peak for the case AV=t shown in Fig. |2] the conduc- 
tance and the phase for this case is presented in more detail in 
a narrow V g range in Fig. |3 We can observe that the partial 
conductances a a and ap are individually small (< 0.25) but 
as the phase between the two channels is they contribute in 
a constructive way, inducing a peak in the total conductance. 

Our results have been obtained solving exactly a relatively 
small cluster and, then, approximately reaching the bulk limit 
using the Dyson equation. It is important to verify that our 
results are robust if the cluster size increases. In Fig. M the 
convergence of the conductance with increasing cluster size at 
AV=0.2t is shown. Remarkably, the diference between L=8 
and L=20 is less than 0.5%, showing that the results presented 
here are nearly exact. The fast convergence is a consequence 
of the local characteristics of the dot, namely all the impor- 
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FIG. 8: Conductance (in units of e 2 /h), phase difference and partial 
conductances for AV=t for the central peak shown in Fig. |2jd). The 
partial contributions to the conductance (a a and op) have the same 
shape as in Fig. |4]but now the phase difference between them is 
zero. This phase difference creates a peak where in the other situation 
appears a dip. 




FIG. 9: Convergence of the conductance with the size of the exactly- 
solved cluster L for the case AV=Q.2t. Dotted lines are for L = 4, 
continuous thin lines for L — 8, and dashed lines for L — 12. Extra 
points for L — 16 and L — 20 are shown with crosses and diamonds, 
respectively. 



tant physics effects are contained in the small cluster solved 
exactly. The fast size convergence shows that even using the 
smallest cluster, the results are already qualitative correct. 

It is important to remark that the cancellations observed 
here do not crucially depend on the couplings U' and J, since 
they are present in our results even for U'=J=0, as shown in 
Fig. [K)] Also for completeness, in Fig. II lithe case of a large 
U'=2t is shown. Even though the situation U'>U is unphysi- 
cal - since the Coulomb repulsion in the same level should be 
larger than for electrons in different levels - it is interesting 
to note that the U' energy plays a role similar to that of Ay. 
Since the energy to have a second electron in (3 is related to 
U' — J + AV, for a large U' (with a J with the same sign as 
considered before) the system reacts as with large AV. Figure 
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FIG. 10: Study of the conductance and total charge when U' and 
J are not active. Shown are results for the case U'=J=0 (solid line) 
compared with the previously shown results at U'=2/3U and J=U — 
U', both at AV=0.35. The effects are qualitative the same in both 
cases. 
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FIG. 11: Conductance (in units of e 2 /h) with large U'=2t. In (a), it 
is shown the case with a ferromagnetic Hund coupling J=\U — U'\ 
and in (b) with antiferromagnetic J=U — U' . Both figures are for 
AV=0. 



^2 contams two differents cases for AT^=0, one of them with 
a ferromagnetic Hund coupling J= \U — U'\>0 (a) and the 
other with an antiferromagnetic coupling J=U — U'<0 (b). 
In the first case, since the Hund coupling is proportional to 
U', we observe the three-peaks structures for the conductance 
with a central-peak associated to the Su=^ states. In the sec- 
ond case, U' and J try to eliminate the <Sd=1 states, and the 
conductance has a similar structure as shown in the previously 
study example with AV=t (see Fig.|2jd)). 



IV. INTEGER-SPIN KONDO EFFECT 

Recently, transport through a dot connected to two leads 
has been experimentally studied as a function of the dot level 
separation (AV), tuned by an external magnetic field. For a 
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FIG. 12: Schematic representation of the different situations that 
could occur in a two-level QD as the gate voltage changes. For more 
details see the text. 



dot with an even number of electrons, it has been shown that 
the Kondo correlations formed between the 5d=1 triplet state 
of the dot and the spins of the leads has a strong influence on 
the current 11 . This problem was theoretically studied using a 
mapping onto a two-impurity Kondo modelA2iii 

Figure [2] shows schematically the different situations that 
could arise as the number of electrons inside the QD is 
changed using the gate voltage. Starting from the upper left 
(a), we have the case where the gate potential V g is above 
the Fermi level and there are no electrons inside the QD. 
Level a is separated from level j3 by the energy AV. De- 
creasing V g , as shown in Fig l^f b). one electron populates 
level a, and the energy to add an extra electron in (3 will be 
AV + U' — J, since an extra energy U' is paid due to the 
Coulomb repulsion and a — J energy due to the spin cou- 
pling. The system now is in the Kondo regime for trans- 
port if T<Tk- Decreasing even further V g , a second elec- 
tron could enter the dot. The spin of this second electron de- 
fines the transport properties since the total spin of the QD 
can be 5d=1 (with a Kondo resonance for temperatures be- 
low ?k(s=i)) or 5*d=0 with no Kondo effect. In the right- 
lower frame (d) the second electron entered in the lower en- 
ergy state a. The energy of this configuration will be of the 
form Es=Q=2V g +U. On the other hand, if the second electron 
enters to the dot with the same spin as the first one forming an 
5d=1 state, we will have the case showed in the left-lower 
frame (c). The energy of this state in the Kondo regime will 
be of the form E s =i=2V g +AV+V-J- T K{s=1 y&. For 
AV>U— £/' + ?k(s=i) we ex P ect a transition from the 5 = 1 
state (right-lower frame) to the 5 = (left-lower frame). 

Note that the spin of the dot 5d is not a well-defined quan- 
tum number, since the QD is not isolated. In the L = 4 
cluster, we found that the ground state used for the calcu- 
lations has a spin 5d=0.7426 (calculated from the relation 
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FIG. 13: Conductance (in units of e 2 /h) for differents J couplings, 
(a) corresponds to the case J=U — U, (b) is at J=0, and (c) is for a 
large negative case J=— 0.5t. The central peak disappears when the 
ferromagnetic spin correlations are eliminated at even larger negative 
J. More detail can be found in the text. 



<5d )ciustcr=5 D (5D + 1)), at AV=0.2 and V g =-0.6 (particle 
hole symmetry). Since the number of particles at this gate po- 
tential is 2, we deduce that the ground state can be expressed 
as a combination of states with 5d=1 and 0, with a weight of 
86 percent for 5d=1- In the case AV=t, we obtain for the 
particle-hole symmetric gate potential a total spin 5d=0.352, 
which corresponds to a 60 percent chance of having 5d=1- 

One way to confirm that the central peak in Fig. |4ja) indeed 
originates in a 5d=1 state is by reducing the value of J, mak- 
ing less probable the ferromagnetic coupling between levels 
a and (3. Figure^]shows the conductance for the case shown 
in Fig. |4ja) plus two differents cases with the same values of 
U, U' and t' but different J's. In (b) with J=0 we observe a 
narrow central peak showing that the effect persists due to the 
Tjc(s=i) energy. For negative values of J, as shown in (c), the 
ferromagnetic state is even less probable and the central peak 
is more narrow. This peak must dissapear for larger negatives 
values of J (i.e. \ J\ >> Tk(S=i))' All this evidence suggests 
that it is the spin-one Kondo effect that is responsible for some 
of the effects discussed here. 



V. UNDERSTANDING THE INTERFERENCE IN THE 
CONDUCTANCE 

The previous analysis showed that the origin of the three 
peaks in the conductance can be traced back to a Kondo effect 
that occurs for 1, 2, and 3 electrons inside the dot. However, 
it remains to be understood what causes the conductance can- 
cellation between them. It is expected that the two paths - 
through a and f3 - are responsible for the interference, but a 
more detailed explanation would be desirable. 
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FIG. 14: Density-of-states of the bare L = 4 cluster, for the cases 
AV=0.05i and V 9 =-1.082i. The Fermi level lies at w=0.0. Note 
the relative change of the weight signs for the poles near u = 0, 
when comparing the Green functions g a N and gprf. 
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FIG. 15: Comparison between the imaginary parts of g aa and g a N, 
for the cases AV=0 and Al/=0.05i, working at V g = -1.082* for 
the L = 4 cluster. The pole in g a N near the Fermi level disappears 
(i.e. carries no weight) at AV=0. 



A. Origin of the it Phase Difference 

To understand how the phase difference arises let us first 
study the cluster Green function. In Fig. ^] the real and 
the imaginary parts of g ar and gp r are shown for the case of 
Ay=0.05t and a gate potential where the relative phase is 7r 
(Vg=0.108i). Since the cluster is by definition a finite sys- 
tem, these two Green functions are simply a collection of a 
large number of poles. For AV=Q, level a is equal to (3 in 
energy, and both Green functions must have poles at exactly 
the same positions. We expect that for a small Ay their fea- 
tures must still remain similar. The most important informa- 
tion arising from Fig. ^3 is that the poles at approximately 
the same energy corresponding to g ar and g@ r near the Fermi 
level (u> — 0) have different signs in their weights. This 
means that those single poles have a phase difference of ir, 
which causes the interference. 

Other special property of the bare (cluster) Green function 
is shown in Fig. ^1 I n this figure, the imaginary parts of 
g aa and g ar for AV=0.Q5t and AV=0 are shown. While for 
AV=0.05i both Green functions have poles at the same posi- 
tions, at AV=0 the weight of the poles near the Fermi level 
for g ar are zero. This is a characteristic of a localized state. 

To understand these two properties we have to analyze the 
non-diagonal Green functions g ar and gp r in their Lehmann 
representation^ 



g ar (iv) = 



(f\c a „\4>l) (<t>l\4M) 



(Ef-Ei) 



iq 



^ (/|C+ (MCaM) | (pj 



[Ey — Ei) 



17] 



gf} r (u) = Y 



(E-p — Ei 



17] 



E 

i 



(/| C + |0 Z ) (frMf) 



(E F - E t 



17] 



(13) 



where \f) and E-p are the ground state and its eigenenergy, 
{|</>i)} is a complete basis of Hamiltonian eigenstates with 
eigenenergies £), and n is an small number (typically 10 -7 ). 

Both Green functions have poles in the same positions 
uj = ±(Ep — Ei). For simplicity, only the first term in Eqs. 
[2Jand^]will be analyzed, but the study presented below is 
qualitatively the same for the second term. As the matrix el- 
ement (4>i\c^ a \f) is the same for both Green functions, it is 
clear that they only can have a phase difference in 7r if the ma- 
trix element (f\c a(J \(f>i) has a different sign than {f\cp a \4>i). 

This change in the sign of the matrix elements can be un- 
derstood based on the "reflection" symmetry between levels 
a and (3. For the case AV=0, the transformation a — > (3 must 
leave the Hamiltonian invariant. The operator r associated 
with this operation commutes with the total Hamiltonian and 
has eigenvalues 1 and —1 (if the states are even or odd under 
this transformation). The complete basis {</>/} used in Eas[T2l 
and in^]can be expanded in two new basis } and 
depending on whether the states are even or odd: 



O r 



(14) 

Or\4>i) = "IO' d5) 

The destruction operators for the a and (3 levels transform as 



0+ C0 a O r 



(16) 



Suppose the ground state is an even state |/ + ) (we confirmed 
this numerically). Then, there are two kinds of matrix ele- 
ments, {f + \c aa \(j>^) and (/ + |c ao -|0j~). These two matrix el- 
ements must be compared with their corresponding matrix el- 
ements, calculated with the operator cp a . Using Easll5land 
I16lwe have 



(17) 
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{f + \c aa \<i>i 



■(f + \<**\ti 



(18) 



Then, some of the poles of Easll2landll3lcarrv weight with 
differents signs. When Ay is different from zero, the weights 
are different but the difference in the signs remains. 



B. Interference between S=l/2 and S=l Kondo states? 

Note that the conductance cancellations mainly occur for 
values of the gate voltage where the charge of the dots is 
changing, for example from 1 to 2. This suggests that an alter- 
native way to visualize the interference process is to imagine 
the ground state in this intermediate regime as a linear com- 
bination of two Kondo states, one corresponding to a spin 1/2 
at the dot and the other to a spin 1 at the same dot. In fact, in 
the important range of gate voltages considered in our study, 
the charge at the dot is not sharply defined, thus this approxi- 
mation is reasonable. The interference could arise from a de- 
structive process between the transport corresponding to each 
of these two nearly degenerate states. This hypothesis will be 
analyzed more extensively in future studies. 



C. Localized States 



Note that the matrix element (f + \c^ a .\4>f) must be zero 
when AV=0, since the operator c+ CT does not change the sym- 
metry of the state \4>j) and the resulting state is orthogonal to 
(/ + 1. This fact produces localized states in the QD due to the 
poles associated with the states \<fi~) (located at uj = Ep — Ei) 
which are present in g aa and gpp (the element (f + \c+ a \(j)f} 
does not appear in these Green functions), but they are not 
present in the Green function that transports electrons from 
one lead to the other (g ar or gi r ). This fact explains the re- 
sults observed in the right frame of Fig. ^] When AV is not 
zero, this is no longer valid since O r does not commute with 

One way to eliminate the symmetry between a and (3 could 
be by using a different diagonal Coulomb repulsion for both 
levels. In Fig. |l6|we can see the conductance, the phase dif- 
ference, and the partial conductances for the case with AV=Q 
and Up = U a + O.K. Note that the first discontinuity remains 
since for one or zero electron in the QD, the different Coulomb 
repulsions U a and Up are not active and the symmetry is not 
broken. In this case all the physics related to the poles associ- 
ated with the symmetry a-f3 discussed above could be applied 
and the discontinuities shown in Fig. |2ja) appears. 
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FIG. 16: Non-symmetric a-/3 case, where the Coulomb repulsions 
in both levels are different i.e. U/3=U a + O.li. Dips in the conduc- 
tance are observed even with A V=0, due to the explicitly broken 
symmetry. More details can be found in the text. 





FIG. 17: (a) Schematic representation of a system of two coupled 
QD, one of them side-connected. In this case, each dot has only one 
level, (b) Equivalent system after the transformation Eqs ll9l and l20l 
The two possible paths in the new basis are shown. 



A. Two Coupled QD 

The first system to analyze consists of two coupled one- 
level quantum dots, one of them side-connected^^, as shown 
schematically in Fig. I17f a). The two dots (A and S) are de- 
scribed by the Hamiltonian H=Ha ts+H\ ca d s +Hi nt , where 



H, 



dots 



H\„t = 



{U/2 nxanxs +V g n\„} 

(7, A— A,S 

£a^C+ CT C sa + h.C, 
a 

*2 ^2 c+ CT [civ + c r<7 ] + h.c. 



VI. CONDUCTANCE CANCELLATION IN SIMILAR 
SYSTEMS 

In this section, we will discuss two other systems where a 
similar interference process is also operative. 



As these are two physically separated quantum dots, we will 
neglect the terms U' and J used in the previous section. Here, 
dot A is the only one connected to the leads, by a matrix ele- 
ment ti. if lead is the same as in previous sections. 

Considering the two dots as a "dimer", i?dots can be rewrit- 
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FIG. 18: Conductance (in units of e 2 /h), phase difference and par- 
tial conductances as a function of the gate potential V g , for the ge- 
ometry shown in Fig^| The parameters used, in units of t (the 
bandwidth is W = 4£), are t a = 0.12, t 2 = 0.25 and U, = 0.5. In 
(a), two zeros in the conductance as a function of V B are shown, (b) 
The phase difference A«I> near these zeros jumps to — n. In (c), the 
partials contributions to the conductance <r_ and <r+ are given. For 
V g — —Ui/2, we can observe a constructive interference process 
where cr_ = <j+ and A$ = 0. 



ten in the basis of bonding and anti -bonding states , 

C+a = -J= [c A a + Cscr] , (19) 
C-a = -J= [c AC r - Cscr] ■ (20) 

With this transformation, it is easy to visualize the two pos- 
sible paths that produce the interference in the conductance 
(shown in Fig. I17> . In the new basis, the matrix element that 
connects the dimer state with the leads is t'=\j \[2 t^- Note 
that the total charge in the two dots does not depend on the ba- 
sis used, since (n Aa + n S(T )=(n_ (T + n +a ). Also note that this 
bonding and anti-bonding states have a diagonal single energy 
of value Vg_=Vg — t a and V g+ =V g + t a . For this reason, the 
matrix element t a plays the same role that 2AV had in the 
case of a quantum dot with two levels analyzed before. 

As in the previous sections, we can define the conductance 
as a function of the partial conductances through channels 
"+" and "— ", as well as the phase difference between them 
as a=a+ + <r_ + ^o + (7_ cos(A$). In Fig. ED we show the 
results for the conductance, phase and partial conductances 
cr+ and cr_ for the case equivalent to that presented in Fig. [2j). 
Here, the dot coupling factor t a used is 0.12<. Again, a three- 
peaks structure in the conductance can be observed. The zeros 
in the conductance appear when the phase jumps to — tt. This 
is a destructive interference process similar to that shown in 
Fig.Hfor AV=Q.2t. 

As in section |m the results presented here correspond to 
the smallest cluster size L = 4. In the final part of this subsec- 
tion the convergence with the size L will be analyzed. In Fig. 
Il~9l the conductance and the charge of the dimer states and 
dots are shown. The charge of the dimer states presents the 




FIG. 19: Conductance (in units of e 2 /h) and charge as a function 
of V g , for the same parameters used in Fig. 1181 Shown are (a) the 
conductance vs. V g ; (b) the charge at the bonding and anti-bonding 
dimer states; (c) the charge at the real dots "A" and "S". 



same behavior as the charge of the two-level dot studied be- 
fore. The charge in the individual dots shows a slow increase 
for dot "A", and a step-function-like behaviour for dot "S". 
This indicates that while the dot coupled to the leads is in a 
transition-valence regime and the charge enters slowly due to 
the Kondo resonance, the other one is in a Coulomb-blockade 
regime. 

Increasing t a the dimer states are splitted and the two-dots 
system can now be represented by one QD with two levels, 
with a large AV~2t a - This is shown in Fig. |20| In the last 
two figures (Figs. I22> the convergence with the cluster 
size is shown for the conductance and the charge. As before, 
a fast convergence is observed. In the case of the charge of dot 
"S", note that the wavy behavior found for the smaller cluster 
is replaced by a "plateaux" when the cluster-size increased. 



B. Two coupled QD with two levels 

In this subsection, it will be analyzed the case of two QD 
in series, each with two levels, as schematically shown in Fig. 
|231 The Hamiltonian that describes this problem includes 



#dots = ^ < U/2niXan iXs + 

i=l,2 [ct,A=q/3 

U'/2 >^ n iaa mp a ' + 

~J C tacr C ioi<7' c tfio' C ^tx + 

<7<7' 

V g [ + (Vg + AV) 

*" [4x*C2X* + cc] , (21) 

cr,\=a/3 
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FIG. 20: Conductance (in units of e 2 /h) for three different values of 
t a - (a) t a =0.12 (reproduced here for comparison), (b) £ a =0.24, and 
(c) £ a =0.48. For t a =0 the side coupled dot is not connected and the 
conductance presents the characteristics of one QD (not shown). In- 
creasing t a , the side coupled dot starts to participate in the transport, 
and interference effects appear. For large t a , the system - repre- 
sented by the dimer states - is equivalent to one dot with two levels, 
separated by an energy 2t a + U. 




FIG. 21: Convergence of the conductance with the cluster size, for 
the case of the side-connected dot. The parameters are the same as 

in Fig, ran 

Hint =f ^2 [ C l\a c l0<? + c 2Xa C r0a + C.C.] , (22) 

where cf Xa creates an electron at the dot i in the level A with 
spin a. The Hamiltonian which describes the leads is the same 
as before. Then, the total Hamiltonian is H^=H dot + -Pleads + 
H lnt . As in the case of two QD with one level each, we can 
rewrite this problem in the dimer-state basis. It is well known 
that in this kind of systems, a splitting in energy between the 
dimer states must occ\x^^^^. This separation is related to 
the t" matrix element. 

As before we present the results corresponding to the small- 
est cluster size. For this case L = 6 that corresponds to two 
sites to describe the "leads" and four levels to describe the 




Vg 

FIG. 22: Convergence study of the charge at dot "S" increasing the 
size of the cluster exactly solved. For the smaller cluster, we found 
important changes around V g =—U/2, but eventually the system con- 
verges to a plateaux for a sufficiently large cluster. The parameters 
are the same as in Fig. 1181 
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3> B 

FIG. 23: Schematic representation of the system of two coupled QD 
with two levels. The two possible paths carrying different phases $ Q 
and <&fj are shown. 



two dots with two levels each. Figure [23] is the equivalent of 
Fig |2] for the case of only one QD. The conductance for sev- 
eral AV's is shown. As in the previous case, at large AV 
there are two splitted structures corresponding to the Kondo 
regime for the 2QD with one level. An intermediate case 
AV=0.9t, where these two structures start to merge, is also 
shown. When AV=Q, similar structures can be observed as 
found in previous cases, where every discontinuity was split- 
ted in two due to the dimer bonding anti-bonding states. 

In Figures [25] and [26] the conductance and the phase for a 
very small Ay=0.01t are presented. In this case the AV=0 
discontinuities are transformed into narrow dips. In Fig. [26] 
the first two dips are shown. They qualitatively correspond 
to the splitting of the first dip in the 1QD case. The phase in 
Fig|26]does not present the abrupt jumps that were found for 
1QD, but the dips nevertheless occur when the phase takes the 
value A$=±7r. 
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FIG. 24: Conductance of two coupled quantum dots, each with two 
levels, for three different cases of AV (large, intermediate, and zero), 
(a) When AV is the largest energy, two structures separated by AV 
appear in the conductance, similarly as for two coupled one-level 
QDs. (b) When AV decreases these two structures start to overlap 
producing dips in the conductance, (c) A more complex conductance 
structure with several discontinuities is obtained in the limit AV=0. 




FIG. 25: Conductance (in units of e 2 /h), phase, and partial con- 
ductances for a small AV=0.01t, in the case of two coupled two- 
levels QDs. Typical dips in the conductance will be analyzed in the 
next figure. Near the gate potential V g =— 0.5, a constructive interfer- 
ence occurs. The phase is zero and the partial conductance becomes 
~0.25, inducing a constructive phase process for electron transport. 



VII. CONCLUSIONS 

In this paper, we have investigated systems where interfer- 
ence effects in the conductance can arise in multilevel quan- 
tum dots. Although there are no obvious different paths, as 
in a AB interferometer, transport using different levels inside 
the QD can produce such interference. Using an exact diago- 
nalization plus embedding technique the phase difference be- 
tween these paths was calculated. When this phase difference 
is ±7r, a destructive interference can appear. An important 
aspect of these results is that the conductance cancellations 



1 - 



C 0.5U 



O 

<3 

-1 
h 4 



" (a) 












; (b) s- 


J] 






V 


V 


- (c) 







-0.1 



Vg 



-0.05 



FIG. 26: First two dips in the conductance for the small AV case 
shown in Fig. 1251 Both partial conductances (a a and 073) have a zero 
near the dips, but not at the same time. Then, the dip in the overall 
conductance is produced by a phase interference process induced by 

A$ = ±7T. 



can occur in the absence of external magnetic fields - usually 
used in AB experiments - and also without introducing at the 
Hamiltonian level any relative phase for the hopping ampli- 
tudes of the states involved. The Coulombic interaction is the 
main responsible for the effect, which can also be visualized 
as an interference between 5d=1 and 1/2 Kondo states. 

Another important novel effect found in this investiga- 
tion is the possibility of having localized states inside the 
QDs. These states do not contribute to the conductance but 
they modify the many-body physics of the systems through 
Coulombic effects. These localized states induce discontinu- 
ities in the conductance vs. gate voltage for particular cou- 
plings, while in the most general case conductance cancella- 
tions are observed. 

The results presented in this paper will hopefully provide 
motivation to experimentalists to refine their measurements to 
search for conductance cancellations. While our conductance 
vs. gate voltage curves should be cautiously used to compare 
against experiments - since the details of the results vary as 
a function of parameters, and in addition most results were 
obtained using only two active levels - the ideas introduced 
here are expected to be robust. When multilevel quantum dots 
with interacting electrons are considered, cancellations are the 
rule more than the exception. 

Experimental realizations of quantum dots with many ac- 
tive levels are possible. Recently, Aikawa et ali 1 observed 
that quantum-dot states can be classified into a small number 
of states that are strongly coupled to the leads, and a large 
number of states weakly coupled to those leads. The mixing 
between these states was found to be of relevance to explain 
transport experiments. It is clear that many theoretical cal- 
culations must be revisited to incorporate the multiplicity of 
levels in dots. This is also important if transport properties 
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through atoms are studied. This is more challenging experi- 
mentally than the study of relatively larger size quantum dots, 
but possible. In the case of atoms, degeneracies are natural 
in open-shell atoms, and a variety of interference processes as 
those observed here are possible. 
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